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I. INTRODUCTION 

The control of light-pulse propagation in matter is a key 
element of optical devices for fundamental science as well 
as for technological applications. In many cases, this is 
made difficult by the presence of competing effects like 
dispersion and absorption. Furthermore, the available 
time for manipulation is often limited by the very high 
propagation speed of light in conventional materials. 

New perspectives in light propagation are opened by the 
observation of long-living coherence effects in optical me- 
dia. The Coherent Population Trapping (CPT) [ll and 
Electromagnetically Induced Transparency (EIT) [2] ef- 
fects have been shown to produce strong modification to 
the properties of optical media. By dressing the mat- 
ter excitations with coherent external fields, a resonant 
probe laser pulse can be made to propagate across an oth- 
erwise strongly absorbing medium at an ultraslow group 
velocity and without being distorted. The incoming light 
is coupled to a dark-state polariton which shows vanish- 
ing absorption and dispersion Q and whose group veloc- 
ity can be controlled via the intensity of the control field 
[4j. This feature has been demonstrated as a tunable 
delay lines for propagating pulses. 

The dynamical modulation in time of the control field 
while the pulse is propagating opens up even richer pos- 
sibilities for light manipulation in the spirit of the so- 
called Dynamic Photonic Structures 5, 6]. For example, 
by completely switching off the control field, the probe 
light can be halted and stored as an atomic (spinlike) 
excitation, and later retrieved after a macroscopic time: 
such light-storage techniques [3, 0] are considered crucial 
tools for all-optical information technologies. A periodic 
dynamical modulation of a spatially homogeneous control 
field can lead to intriguing phenomena such as frequency 
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triggering in time of the EIT band Q . By adiabatically 
raising the resonance frequency, it is also possible to shift 
the energy of a propagating polariton Q . A nonadiabatic 
variation of the control field has been proposed as a tool 
to compensate the pulse broadening at the exit of a de- 
lay line [13 or after retrieval of a previously stored light 
wave packet [ll[ . Extremely fast modulations of the con- 
trol field have been anticipated to produce a substantial 
dynamical Casimir emission (l2|. 

Spatially modulated EIT media have been proposed and 
used for light-stopping applications, via the creation of 
tunable stop bands within the EIT window [H, [TH . Mu- 
tual interaction of several moving spin coherence gratings 
has been demonstrated as an efficient way to stop two- 
color light and to perform wavelength conversion [14j . 

In the present article we investigate the effect on ultra- 
slow light propagation of a combined spatial and tempo- 
ral modulation of the EIT medium. A pioneering study 
in this direction was reported in flH HH ] where the simul- 
taneous propagation of both control and probe pulses is 
investigated: a ramp of the control field in an otherwise 
homogeneous medium induces different propagation ve- 
locities in the different parts of the probe pulse. This 
is predicted to result in a controllable reshaping of the 
probe profile. 

Dilute ultracold gases are among the most promising me- 
dia for EIT applications. Both slow light and light stor- 
age have been experimentally realized in this systems 
|17| . As compared to condensed matter systems, atomic 
gases have the advantage of showing extremely narrow 
linewidths, which can be exploited to obtain ultraslow 
group velocities: Doppler and collisional dephasing pro- 
cesses which destroy the coherence are in fact strongly 
reduced in cold samples. Furthermore, the achievement 
of Mott-insulator (MI) phase [l8| suggests the possibility 
of obtaining even longer coherence times thanks to the 
gap in the many-body excitation spectrum: the first ex- 
perimental realization of EIT in a MI has been recently 
reported for light-storage purposes [l9j . 
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Most of the theoretical works on slow light propagation 
in atomic samples were focused on the simplest case of a 
homogeneous atomic medium with some boundary con- 
dition. Unfortunately, the typical size of atomic samples 
is often small as compared to the duration and waist of 
the probe pulse, which calls for a more complete theoret- 
ical description. In this article we present a model that 
is able to include the spatial inhomogeneity of a system 
and therefore to describe the propagation dynamics at 
the interface between vacuum and the EIT medium. 

In particular, we show how one can take advantage of 
the interfaces to manipulate the wave-packet shape by 
means of a dynamical modulation of the control field in- 
tensity. The multilayer structure offers in fact the possi- 
bility of spatially engineering the radiation-matter inter- 
action. While the usual schemes based on a single control 
field pulse in a homogeneous system [l5[ limit the manip- 
ulation to a single interface, our proposal can explore a 
much wider range of configurations: increasing the num- 
ber of interfaces allows a variety of modulation protocols 
to be developed, leading, for example, to the switching 
of a single pulse into a train of separated signals. Fur- 
thermore, the reduced optical depth of each layer may al- 
low for a more efficient modulation Q . The short pulses 
that can be tailored with our proposed technique may 
be of interest in view of creating polariton states with a 
Dirac-like dispersion [2(| ■ Even though our theoretical 
description is based on semiclassical wave equations, our 
conclusions directly apply to the quantum processing of 
single photon wavefunctions and therefore can have an 
importance for optical quantum computing applications. 

The article is organized as follows. In Sec. II we re- 
view the Maxwell-Bloch formalism for describing light 
propagation in an atomic medium and we introduce the 
modified Slowly Varying Envelope Approximation to in- 
clude the spatial inhomogeneities. The standard polari- 
ton picture for slow light propagation in homogeneous 
EIT media is reviewed in Sec. III. These concepts are 
then used in Sec. IV to investigate the simplest modu- 
lation schemes based on the idea of EIT chain. In Sec. 
V, we derive an effective equation for the pulse propaga- 
tion. This simplified formalism is then used in Sec. VI 
to propose and characterize some specific protocols to be 
implemented in cold atom systems. In Sec. VII, we draw 
conclusions and we sketch the perspectives of the work. 



II. THE FORMALISM 

This article is devoted to the theoretical analysis of the 
propagation of a light pulse in a time-dependent and 
spatially inhomogeneous medium with sharp interfaces. 
Furthermore, the pulse carrier is fixed in a region where 
the optical response of the medium is strongly dispersive. 
This requires that the formalism explicitly includes both 
the atomic degrees of freedom in the presence of a time- 
dependent dressing field, and the spatiotemporal dynam- 



ics of the propagating pulse. In this section we introduce 
a modified version of the well-known Slowly Varying En- 
velope Approximation which is able to address all these 
features in a single unified treatment. 



A. Optical Bloch equations 

We assume the internal dynamics of the atoms to be re- 
stricted to three levels in the A configuration sketched in 
FigQJa): the ground state \g), a long-living metastable 
state |m), radiatively decoupled from \g} 7 and an op- 
tically active excited state |e), with a lifetime 7" 1 . 
The lifetime of the \m) state is assumed to be much 
longer than any other time scale in the system: 7 m * 3> 
(7(T 1 ) tmax) where t max is the time of observation of the 
system. The atomic ensemble is assumed to be initially 
in the ground state. 

Two lasers are shone on the atoms: the control (or dress- 
ing) field is a strong laser beam coupled to the \m) — > \e) 
transition of carrier frequency uj c and (time-dependent) 
Rabi frequency 51 c . This field is considered to be an ex- 
ternal field that can be controlled at will. The focus of the 
present work is centered on the propagation of another 
laser pulse of carrier frequency uj p close to resonance with 
the \g) — ¥ |e) transition of dipole moment d eg , and has 
an electric field 

E{x,t) = £ (x,t) e'^r* +c.c. (1) 

This probe field is assumed to be weak enough to be 
within the linear regime where a linearized version of the 
optical Bloch equations 0, Hl| can be used. Within the 
well-known rotating-wave approximation, only quasireso- 
nant terms for the atom-laser interaction are considered. 

In order to eliminate the carrier frequency, we can rewrite 
the optical Bloch equations in terms of the reduced co- 
herences p eg = p eg e""»* and p mg = p mg e l (" P - w c)t. 

dpeg (le . \ . . d ge £ ,O c . 

-gf = - \y + e ) peg + l ~h — t Y Pmgi ^ ' 

dPmg /7m . . r \ - . - h 

—gf = - {— + WflJ Pmg - J yfe- ( 2b ) 

Here, S e = uj e — uj 9 — oj p is the one-photon detuning of 
the probe field from the \g) — > \e) transition, while Sr = 
oj m +L) c — LJg — cjp is the detuning of the probe and control 
beams from the so called Raman two-photon transition 
connecting the ground to metastable states. 

As the transition \g) — ¥ \m) has a vanishing dipole mo- 
ment and the coherence p ern is of higher order in £, the 
atomic polarization can be written in terms of the coher- 
ence Peg only: 

P(x,t) = n(x)d eg p eg (x, t) + c.c. (3) 

Here, n(x) is the atomic density. 



3 




FIG. 1: (a) Scheme of the three-level A configuration of atomic levels and laser beams, (b) and (c) Frequency dispersion 
and decay rate of the three polariton branches (green, upper; blue, dark; red, lower). Parameters of the system: \f~D = 0.1, 
fi c /u; e g = 0.04, 7 e /u eg = 0.01, resonant dressing ui c — uj e —U! m . The solid lines are the MSVEA predictions using uj p = w eg and 
and assuming an Erf shape for the f(k) function with a half-bandwidth of 0.5 ui eg . The dashed line is the exact EIT dispersion. 



B. Modified slowly varying envelope 
approximation (MSVEA) 

In the present work we restrict our attention to a one- 
dimensional geometry at normal incidence. As the dif- 
ferent polarizations of electromagnetic (e.m.) field are in 
this case decoupled, the vector nature of Maxwell equa- 
tions disappears and one is left with a scalar problem for 
each component: 



dx 2 



c 2 dt 2 



E(x,t) 



d 2 



(4) 



Here P(x, t) is the polarization Q of the atomic medium. 
The constants c and no are, respectively, the velocity of 
light in vacuum and its magnetic permeability. 

We assume that the pulse envelope £{x,t) varies on 
a characteristic time scale much slower than the car- 
rier frequency u) v . Under this approximation, we can 
perform a modified slowly varying envelope approxi- 
mation (MSVEA) and neglect the second-order time 
derivatives of the envelope. Assuming the same form 
P(x, t) — V(x, t)e~ tulpt + c.c. for the atomic polarization, 
the Maxwell equation Q can then be written in the form 



d*_ 

dx 2 



2 1 



£ = -poWpP. 



(5) 



Here we have neglected the first-order derivative in time 
of the polarization envelope which is proportional to 
second-order time derivative of the electric field 1221. 



It is important to note that, differently from the conven- 
tional SVEA discussed in textbooks (e.g., [23j]), we have 
not separated out the spatial part of the envelope from 
the carrier and we have retained all the derivatives of the 
field with respect to the spatial coordinates. This crucial 
feature of the MSVEA equation ([5]) is required when one 
investigates configurations involving abrupt jumps in the 
spatial distribution of atoms n(x). 

The optical polarization Q arising from the optical 
Bloch equations has to be plugged back into the Maxwell 
equation After choosing an appropriate normaliza- 
tion, the resulting Maxwell-Bloch (MB) equations can be 
cast in the form 



d£ 
dt 

dt 
dt 



(6a) 



(y + P<^9 + i^f~D£ - i-j-Pmg, (6b) 



In particular, ui p is used as the unit for frequency, and 
k p = {ojp/c) for the wave vector. The electric field is 
measured in terms of £o = y / 'nhuip/2eo. The physical 
meaning of this choice is related to the energy density in 
the system: the energy density associated to the atoms 
is Wat — nhu! p \p eg \ 2 , while the energy in the e.m. field 
is [24j W e , m , = 2€ £q\£\ 2 : £ is then the electric field as- 
sociated to a polariton in which the excitation is exactly 
shared between atoms in the excited state and photons, 
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\£f = \p eg f = l,W a t = W e . m .. 

The strength of the light- matter coupling is quantified by 
the adimensional parameter 



medium are respectively given by the dispersion laws 



D = 



(7) 



In the case of a homogeneous system of N atoms, the D 
coefficient is straightforwardly related to the atom-field 
coupling constant g as defined in Q, 



DloI = g 2 N. 



Features and limitations of MSVEA 



(8) 



Before proceeding, it is important to assess the features 
and limitations of the MSVEA approach that we intro- 
duced in the previous section. This approximation leads, 
in fact, to Eqs. ([3]) and © that differ from the standard 
formalism used for EIT-related problems and offer im- 
portant advantages for the specific problems under con- 
sideration here. 

In the absence of atoms the MSVEA equation ([5]) with 
V = gives the following approximate dispersion for the 
free e.m. field: 



;(fc) = 



2u; K 



(9) 



On one hand, this dispersion is able to simultaneously 
describe both the forward (k > 0) and the backward 
(k < 0) propagating photons. This will be useful to han- 
dle reflectivity problems without the need for a coupled 
mode theory. In particular, this approximate dispersion 
correctly reproduces the exact one in the neighborhood 
of k = ±k p , for example 



u)(±k p ) 
duj 
dk 



c\k p \, 
±c. 



(10a) 
(10b) 



On the other hand, the deviation from the linear disper- 
sion of light that is due to the curvature of the approx- 
imate dispersion (|9]) is responsible for a spurious wave- 
packet broadening. However, this effect start to be im- 
portant over propagation lengths that are much longer 
than the ones under investigation here. 

At the interface with a generic semi-infinite medium of 
linear susceptibility x(w), the reflectivity of a monochro- 
matic wave at normal incidence can be straightforwardly 
calculated from ([5]) as 



i - k'/k 



(ii) 



l + k'/k 

where the MSVEA wave vectors in vacuum and in the 



ck = uip J 1 + 2{u) — U! p ) /ujp, 



(12) 



ck' = oj p \Jl + xM + 2(w - ujp)/LOp. (13) 

Provided the frequency u is close to the carrier ui p , the 
approximate reflectivity (|11[) is accurate up to corrections 
of the order (u— aj p )/ujp. This condition is well satisfied in 
an EIT medium in the frequency region around resonance 
as the light propagation is dominated by the frequency 
dispersion of the susceptibility x(w). 

To conclude this section, it is important to note that at 
the level of the MSVEA approximation one is allowed to 
replace the c 2 k 2 term in ([9]) with a generic function f(k) 
that satisfies the conditions (fT0|) . This feature is of great 
interest when one is to numerically solve the set of Eqs. 
([6]). In fact, a proper choice of f(k) makes it possible to 
significantly increase the time step of the simulation [25[ : 
in the remainder of the article, we will use the Erf-shaped 
form 



/(*) = <4 



1+Erf 



(14) 



which results in a bandwidth of the order of uj pi wide 
enough to avoid introducing spurious physics in the fre- 
quency region of interest. The resulting dispersion for 
the EIT system is shown in FigJU which is in excellent 
agreement with the exact one in the frequency region of 
interest. In particular, the choice of a linear f(k) in the 
vicinity of uj p allows to suppress the effect of the spurious 
dispersion of the wave packet that would be otherwise in- 
troduced by the mSVEA. We have also checked that the 
numerical results presented in the remainder of the arti- 
cle do not depend on the specific choice of f(k). 



III. POLARITONS IN A HOMOGENEOUS AND 
STATIC EIT MEDIUM 



Before entering in the discussion of the complex features 
of the propagation in spatially inhomogeneous and time- 
dependent media, it is useful to briefly review the main 
features of light propagation in the simplest case of ho- 
mogeneous and static EIT medium. 

At the level of linear optics considered in the present ar- 
ticle, the response of a stationary medium is summarized 
in its susceptibility x(w). In the case of a EIT medium, 
this is straightforwardly calculated from the Bloch equa- 
tions [(16bl) and (l6cl)l to be 



X(w p ) = 2Duj p 



x ■> 
o e . - i — 



Sr - i(7m/2) 



(15) 



Once plugged into the Maxwell equation, this form of 
x(w) gives the three polariton branches [H, [l2| that are 
shown in Fig. [lib). I n this figure, as well as in the rest 
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of the article we focus our attention on the resonant case 
uj c = uj e — cu m . The width of the Rabi splitting between 
the polariton bands at resonance is fixed by the light- 
matter coupling to (2\Du)p). Note how the Erf-shaped 
form of f(k) chose for the calculation does not affect the 
polariton dispersion in the region of interest. 

Throughout this article, we are mostly interested in the 
central band, the so-called middle polariton or dark po- 
lariton (DP) Q. The width of this DP branch is pro- 
portional to the dressing Rabi frequency Q c . Exactly on 
Raman resonance 6 e — 0, the group velocity is given by 



dx 



2 duj n 



1 + 



(16) 



(n c /2) 2 



In contrast to standard optical media where an optical 
resonance is generally associated to a high interface re- 
flectivity and a strong absorption [t| [26H-28J , three- level 
systems under a coherent dressing allow for a very slow 
group velocity, a negligible interface reflectivity, and neg- 
ligible absorption Q. This striking effect is due to the 
vanishing value of the susceptibility x( w ) exactly on res- 
onance combined with a strong value of its frequency dis- 
persion. Spontaneous emission from the excited state is 
suppressed thanks to the destructive interference between 
the different excitation paths, while dephasing between 
the ground and metastable states is generally very small 
in cold atom systems, 7 m ~ 0. 

For a polariton wave at a wave vector k, the group ve- 
locity and the lifetime are in fact related to the relative 
weights of the radiation and matter excitation compo- 
nents that are obtained by diagonalizing jj| the set of 
Eqs. @: 



V9r(k) ~ c Wf 



\£{k)Y 



7 (fc) 



7e 



\Peg{k)\ +\Pmg{k)\ 
\Peg(k)\ 2 

\£(k)\ 2 + \p eg (k)\ 2 + \p mg (k)\ 2 ' 



(17) 
(18) 



The prediction (fTBj) for the decay rate of the different 
bands is plotted in FigQJc). As expected, the decay rate 
is exactly vanishing on resonance (for 7 m = 0) and grows 
quadratically in the wave vector k: 



j(k) = -Im 



' d 2 L0(k) 




dk 2 


S R =0_ 



{k kp) 



(19) 



To leading order, the second-order derivative of the EIT 
dispersion is then purely imaginary f22j ]: 



a u 
dk 2 



c <20) 

For a sample length L, this recovers the well-known @, Q 
Gaussian transmittivity window of width 



Aojtr 



1 (Oc/2) 2 

^ y/le DU 2 



(21) 



IV. MULTILAYER SYSTEM: THE EIT CHAIN 

We consider a pulse of light launched into a layered ge- 
ometry consisting of several atomic EIT media separated 
by empty regions of space. A pictorial view of the layered 
system under consideration -the EIT chain- is shown in 
Fig 12] We assume the atoms in the different EIT layers 
to have the same Raman frequency. The probe pulse car- 
rier is taken exactly on Raman resonance and the pulse 
bandwidth is assumed to fit within the EIT frequency 
window. The propagation of the pulse across the system 
can be simulated using the MB formalism © described 
in the previous sections with a spatially dependent D(x) 
and Q c (x). 

Two cases can be distinguished: a static case where all 
the parameters describing radiation-matter interaction 
remain costant in time, and a dynamic configuration 
when some parameter is varied in time while the pulse is 
propagating across the system. In what follows, we will 
concentrate our attention on a time modulation of the 
control field amplitude f2 c . Other dynamic schemes may 
involve a modulation of the atomic resonance frequency 
and/or of the dressing field frequency @, @|. 






Probe 



Control 



FIG. 2: Pictorial scheme of a double-layer EIT chain. 
The shown configuration with orthogonal probe and control 
beams ,17] allows for independent modulation of fi c in the 
different layers, but the predictions of the present article are 
straightforwardly extended to other geometries. 



A. Static case 

In a static situation, the pulse propagates across the 
whole multilayer structure with negligible distortion. 
The propagation time is equal to the sum of the thickness 
of each layer divided by the corresponding group velocity. 

While the spatial shape of the emerging pulse remains 
unchanged with respect to the incident pulse, the shape 
inside the structure is significantly modified. As shown 
in FigEJa) , a pulse entering an EIT layer is in fact spa- 
tially compressed by v gr /c as a consequence of the re- 
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duced group velocity and a reversed process takes place 
when the pulse leaves the layer, which restores the ini- 
tial shape. Correspondingly, the envelope of the electric 
field remains continuous at all interfaces, but its deriva- 
tive has discontinuities proportional to the group velocity 
mismatch between neighboring layers. 

An important distinction has therefore to be carefully 
made in the notation: the length of the wave packet 
within the EIT layer will be denoted by a x = at v gr while 
its length in vacuum will be denoted by a x — at c. 



tional to the modulation rate of the control field ampli- 
tude d£l c / dt. This matrix element is to be compared to 
the frequency splitting between the bands, which leads to 
the following quantitative criterion for adiabaticity [29i — 

a 

< 22 » 

2. Interface 



B. Dynamic case 

The easy tunability of the properties of the dressing 
field together with the slow propagation of the DP al- 
low for an efficient dynamic modulation of the propagat- 
ing pulse. Taking advantage of spatial inhomogeneities, 
the dynamic EIT chain can be the paradigm for a new 
class of dynamic photonic structures @ which are able 
to perform a simultaneous spatial and temporal modu- 
lation of the wave-packet profile. The mechanism of the 
pulse modulation in complex geometries can be under- 
stood in terms of two basic building blocks: the homo- 
geneous layer, when the pulse is completely contained in 
a single homogeneous EIT layer during the whole modu- 
lation sequence, and the interface, when the modulation 
takes place while the pulse is overlapping two neighboring 
layers. 



1. Homogeneous layer 

When the whole pulse fits into a homogenous EIT layer, 
its dynamics can be easily understood within the polari- 
ton picture discussed above and it is fully determined by 
the conservation of wave vector k. 

In the limit of a very slow modulation, the pulse adiabat- 
ically follows the evolution of the DP state Q and propa- 
gates at the instantaneous group velocity with no change 
in its length. The peak electric field follows the magni- 
tude of the photonic component (|17|) and is determined 
by the ratio of the final to initial group velocity /v l . 
In particular, it does not depend on the temporal shape 
of the modulation. Examples of such modulations are 
illustrated in Fig. [3jb). In the case in which the group 
velocity modulation is brought back to the initial value 
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FIG. 3: Examples of modulation of a pulse via a space- and 
time-dependent EIT. The propagation of the pulse is calcu- 
lated by solving the MB equations ((6j). In all of the panels, 
the pulse moves from left to right and the shaded regions cor- 
respond to a EIT medium. The dashed curves show the initial 
pulses, (a) Propagation across a static interface. Solid line, 
pulse shape while being spatially compressed in entering into 
a medium with v gr — 0.11 c; Dotted line, pulse shape once 
it has completely entered the EIT medium, (b) Wave-packet 
propagation in a homogeneous EIT layer. Dotted line, prop- 
agation without modulation; solid line, result of a slow down 
ramp with v gr = v l r /2; dot-dashed line, result of a speed-up 



v l the pulse emerges with an unchanged profile ram P v lr = 1-8 v % gr . (c) Exit from a EIT medium into vacuum 



and the layer can be considered as a very compact, yet 
programmable delay line. It is interesting to note that 
the dynamical modulation of the polariton branches en- 
sures that, if the pulse fits into the EIT window at the 
entrance, then it will fit into it at all later times [H, Q . 

The main effect of a finite ramp time r is to couple the DP 
to the upper and lower polariton branches at the same k. 
The effective matrix element of this coupling is propor- 



Dotted line, propagation without modulation; solid line, re- 
sult of a slow-down ramp with v gr = v gr /2. All the panels 
have been calculated for a pulse carrier exactly on Raman 
resonance S e = 0, a ramp time uj p t — 100, and a pulse width 
Lo p (j t — 400. Material parameters: D = 0.01, 7 e = lO^ujp. 



A dynamical modulation taking place while the pulse 
overlaps an interface provides a simplest way of reshaping 
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the pulse: only the part of the pulse which is located in- 
side the EIT layer is in fact affected by the modulation of 
the dressing field. In contrast to the the spatially homo- 
geneous case considered earlier in this article, the shape 
of the emerging pulse now strongly depends on the details 
of the modulation ramp even in the adiabatic limit. This 
crucial fact is illustrated in FigJ^c): the group velocity 
of a EIT medium is reduced while a pulse is exiting into 
vacuum. 

The results can be understood by isolating three portions 
of the pulse: the first part is already in vacuum when the 
modulation begins, while the last part is still in the EIT 
medium when the modulation is completed. 

The first part is therefore not affected by the modula- 
tion, while the electric field amplitude of the third part 
is homogeneously lowered. When also this part of the 
pulse eventually exits into vacuum, its spatial length is 
stretched out by a factor c/v]j r . As compared to the front 
of the pulse, the stretching is larger by a factor v l gr /v^ r , 
which results into the strongly asymmetric pulse shape 
that is visible in Fig|3]Jc). 

Finally, the modification of the central part of the pulse 
depends in a nontrivial way on the details of the ramp. 
For a fast [but still adiabatic as compared to the inter- 
band splitting, according to ([22]) ] modulation, the first 
and third parts of the pulse are connected by a sharp 
jump in the electric field amplitude. For slower ramps, 
this jump is replaced by a smoother crossover. 



3. Defect 

The physics of a defect geometry can be understood along 
these same lines. Two regions of a homogeneous medium 
are separated by a thin layer with a different group veloc- 
ity. The thickness Ld of the defect region is taken to be 
small as compared to the effective length of the pulse in 
this layer. For the sake of simplicity, we restrict our at- 
tention to the simplest situations of a EIT medium with 
a vacuum defect and of a single EIT slab in vacuum. 

The case of a vacuum defect is illustrated in FigfJ] In 
particular, we consider the long pulse limit a x 3> Ld in 
which the pulse amplitude can be considered as almost 
homogeneous across the defect. The temporal shape of 
the modulation of Q c is shown in the inset: its time scale 
is assumed to be fast as compared to the pulse duration, 
but still slow as compared to the interband adiabaticity 
requirement (I22[) . Vg r are the maximum and minimum 
values of the group velocity; t s is the interval between 
the two ramps, that is, the storage time. 

The modulation takes place in three stages: the slow- 
down ramp, the storage time and the speed-up ramp. 
When f2 c is first decreased, the pulse intensity in the 
atomic medium is correspondingly reduced by a factor 
v gr/ v gri w hile the amplitude of the central part remains 
unchanged as it is sitting in vacuum [FigHJa)]. When 
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FIG. 4: Modulation of a wave packet using a vacuum de- 
fect. In the EIT medium, the group velocity is decreased from 
Vg T = 0.11 c to v~ r = 0.02 c and then increased back to Vg^ r as 
shown in the inset. (Main panel) Pulse at the beginning (dot- 
ted blue line) and at the end of the process (solid blue line); 
comparison with the results of the continuity equation for po- 
lariton flow (|24[) (black dot-dashed line) and the numerical so- 
lution of the effective equation (I29p (red dashed line). (Inset) 
Temporal dependence of Q c . (Bottom panels) Three snap- 
shots during the propagation time, (a) Slow-down ramp; (b) 
storage time; (c) speed-up ramp. Parameters: j e — 10 -3 Wp, 



time uipT = 100, storage time ui p t s = 60000. 



6400, ramp 



this part of the pulse re-enters the EIT medium, it re- 
sults in a spatially compressed, narrower peak of width 
LdV~ r /c [FigQJb)]. The part of the pulse that crosses 
the defect during the storage time does not experience 
any distortion: in FigJ^b) this part lies just behind the 
narrow peak and is T s v~ r — Ld v~g r /c long. The final speed- 
up ramp which restores the group velocity to its initial 
value v~^ r is responsible for an increase of the electric field 
amplitude in the EIT medium. This eventually results 
in a hole being left imprinted in the pulse profile corre- 
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spondingly to the vacuum layer [Fig|4jc)]. Once it has 
re-entered into the EIT medium, the length of this hole 
is equal to LdV+ r /c. 

This vacuum defect geometry is then the simplest exam- 
ple of a nontrivial modulation of the wave-packet profile: 
note that the resulting profile is very different from the 
one obtained in a homogeneous geometry, where the sec- 
ond ramp would simply compensate the first one. 
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FIG. 5: Snapshots of the pulse profile during a light-storage 
process in a single EIT layer. fi c is modulated in time from 
Q t = 0.07 Ld p to £7(7 = with the same shape as shown in 
the inset of FigUJ The storage time is ui p T a — 1350. The 
initial group velocity in the EIT layer is v gr = 0.11 c. The 
pulse has a Gaussian shape with k p a x = 540 in vacuum and 
the defect has a length k v Ld = 10. (a) Initial pulse (dashed 
line), pulse hitting the defect and splitting immediately after 
the stopping ramp (solid line), and counterpropagating wave 
packets during the storage time (dotted line), (b) Emerging 
wave packets after the retrieval ramp. The arrows indicate 
the propagation directions of each pulse. Blue lines in panels 
(a,b) are in the absence of spontaneous emission 7 e = 0. The 
green line in panel (b) is for f e = 0.07 oj p . 



The opposite case of an EIT layer in vacuum is presently 
of great experimental interest for light-storage purposes 
[3, |3] • The idea is very simple: by switching off the 
control field fl c while the wave packet is inside the EIT 
medium, the DP is fully mapped into a metastable co- 
herence p mg . As this has a vanishing group velocity, it 
can remain stored in the atoms for macroscopically long 
times. When the control field ft c is switched on again, 
the wave packet is retrieved. The main limitations to the 
efficiency of a light-storage process originate from sponta- 
neous emission processes from the excited state (a finite 
p eg component is always present for any pulse of finite 
duration) , leakages due to the finite optical depth of sam- 
ple as the usually considered systems are shorter than the 
effective length of the pulse, and ground-state decoher- 
ence j m > [!, 0] ■ Even though this last effect sets the 



ultimate limit to the performances of light-storage exper- 
iments, for the parameters considered in the present work 
its effect is negligible as compared to the other processes. 

A situation similar to the experimental realization in (l9j 
is simulated in Fig[5] because of the limitations in the 
numerical solution of the MB equations, we have been 
forced to consider a EIT medium with a much bigger 
Vgr/c. Anyway, a suitable rescaling of all other parame- 
ters makes it possible to observe the same dynamics. 

During the storage time, fi c is brought to zero with the 
same temporal profile as shown in the inset in Figj4] 
again, this process consists of three stages. During the 
stopping ramp [FigJSJa)], the signal is cut into three 
parts. The front part, which has already crossed the de- 
fect, is not affected by the modulation and keeps on prop- 
agating almost undisturbed. The central slice that was 
contained in the EIT layer at the stopping time remains 
coherently stored in the medium as an atomic polariza- 
tion. The back part of the pulse hits the medium when 
this is no longer optically dressed and its fate strongly 
depends on the spontaneous emission rate from the |e) 
state. 

If spontaneous emission is negligible, this part of the 
pulse is reflected back albeit in a strongly distorted 
way 0, this is for example, the case of a Mott- 

insulator state of two-level atoms where spontaneous 
emission is suppressed as a consequence of the ordered 
lattice structure 



32] 



and information on the interface 
structure can be inferred from the shape of the reflected 
wave packet. On the other hand, as shown by the green 
(light gray) line in FigJSJb), the reflected wave packet is 
almost completely absorbed in the presence of a signifi- 
cant spontaneous emission rate. 

When the retrieval ramp is finally applied, the excitations 
stored in the EIT layer are free to propagate out of the 
atomic medium. The efficiency of this retrieval process 
(defined as the ratio of the intensity of the retrieved pulse 
to the intensity of the incident pulse) is around 15% for 
this simulation, which qualitatively agrees with the esti- 
mation in (l9j . As shown in FigJSJb), this retrieved wave 
packet is weakly affected by spontaneous emission. 



POLARITON FLOW: EFFECTIVE 
DESCRIPTION 



The MB formalism © gives a complete picture of the 
pulse propagation which is able to take into account in- 
terband transitions as well as reflection at the interfaces. 
As the solution of the three coupled equations is time- 
and memory-consuming, it quickly becomes unfeasible 
for growing values of the velocity mismatch between the 
different media. 

For this reason, an effective approach able to investigate 
the ultraslow light regime can be of great interest. As 
the dynamic modulation of the pulse does not result into 
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an increased reflection as compared to the static case, 
the propagation dynamics can be described in terms of a 
single continuity equation for the polariton flow. In par- 
ticular, both absorption and the spatial inhomogeneity of 
the structure are fully included in this model equation. 



A. Continuity equation 

In terms of the polariton density n p (x,t), the continuity 
equation describing conservation of the total number of 
polaritons reads 



d 9 . 

—n p [x,t) + — (n p (x,t)v gr (x,t)) 



0. 



(23) 



In the static situation (v gr (x,t) — Vg r (x)) 7 also the local 
polariton flux n p (x, t)v gr (x) remains constant. It is use- 
ful to rewrite this equation in terms of the electric field 
intensity / = \6\ 2 corresponding to the polariton flux 
instead of the polariton density. For static multilayer ge- 
ometries which are carachterized by abrupt changes in 
the polariton velocity, the polariton density shows, in 
fact, discontinuities while the electric field remains con- 
tinuous, even at the interfaces. Taking into account the 
photonic weight and the group velocity (fT7|) . Eq. ([23|) then 
becomes 



m 

dt 



v 



I dv 



gr 



dt 



(24) 



The left-hand side of the equation contains the propa- 
gation terms for the static situation. The general so- 
lution is a mixed translation and dilation of the start- 
ing pulse I°(x) according to the trajectories in space- 
time which are the solution of the Cauchy problem: 
{£ — v gr (t;),£,(t) = x) [HHU. The specific solution clearly 
depends on the geometry. The right-hand side of (f2"4")l is 
instead responsible for the amplitude variation in the dy- 
namic case. 

As a first example, this continuity equation model can be 
applied to the defect geometry of FigJH 

Vgr(x,t) = v(t) {6{-x) + 9{x ~ L d )) + c9{x)6{L d - x), 

(25) 

where the beginning of the defect is located at Xd = 0. 
For times longer than t — (Ld/c), the solution for the 
electric field intensity is 



2° (a 
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Here t d (x) and to(x) are the instants of time at which 
the point of the wave packet which is located at x at the 
time t has passed through x — L d and x — 0, respectively. 
They can be found by the conditions 



L d 
--t d - 



■A 

hi 



L d /c. 



where we have defined 



l b a = / v(t')dt'. 

J a 



(27a) 
(27b) 



(28) 



This analytic solution gives the black dash-dotted curve 
in the main panel of Fig01 



B. Effect of losses 

Even if the carrier frequency u> p sits exactly on Ra- 
man resonance, the finite time duration of the wave 
packet requires including absorption for the tails of the 
wave- vector spectrum [33J. This leads to a finite and 
momentum-dependent decay rate for the polaritons ac- 
cording to (|18p . Taking inspiration from the approxi- 
mated form (fin]) of the decay rate, a simple diffusion 
term can be used to model losses (34|. The propagation 
Eq. for the intensity / then becomes 



~dt 



dl_ 

dx 



U gr dt 



dx dx' 



where 



Jgr 



C7e 



(30) 



is the diffusion coefficient. As shown in the main panel 
of Fig. HI this effective model (red dashed line) cap- 
tures all the features of the full MB calculation (blue 
solid line) with a good accuracy [35]. In particular, the 
effective model is able to perfectly reproduce the height 
of the peak that was instead overestimated by the simple 
diffusionless continuity equation (|24l) (black dash-dotted 
line) . 



VI. MANIPULATION SCHEMES IN COLD 
GASES 



The EIT chain introduced and discussed in the previous 
section can be implemented experimentally using clouds 
of ultracold atoms as the EIT media. Optical fibers can 
be used to fix the optical distance between the EIT layers 
36]. Other possible material realizations involve solid- 
f s&c materials |2|. 

Realistic values for the system parameters can be ob- 
tained from [T3, H3- As a typical example, we con- 

19 , 



(26) sider a cloud of Na atoms of density n = 8 x 10 m 
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For the optical transition, we use the D2 line. As the 
ground state, we use the \g) — \3Si/2,F = I, Trip = — 1) 
sublevel; as the metastable state, we can use \m) = 
|3S'i/2,-F 1 = 2,mp — —2), and as the excited state we 
can use |e) = \3P 3 / 2 ,F = 2, nip = —2). In this case, 
we have u eg = (2tt)508 THz, d eg = 1.5 x 10" 29 Cm, 
7e = (2tt)10 MHz and D = 3 x 10~ 9 . For a control 
field of Rabi frequency ft c — (27r)17 MHz, a group ve- 
locity v gr — 10 _7 c is obtained along with an absorption 
coefficient T) = 6 x 10 -7 ijj eg jk\ g . The parameters of the 
recent experiment are not much different. From here 
forward, physical units are used in the figures. 

We have used the effective Eq. ([21))) to simulate several 
simple geometries. We first address the single-layer case 
and then we use the results to investigate more compli- 
cated multilayer structures. 



A. Single layer 

We consider a single EIT layer where we inject a Gaussian 
pulse. As discussed in Sec. IV, the effect of the modula- 
tion depends on several parameters: a crucial quantity is 
the ratio R = L/a x between the layer thickness and the 
pulse length. Here, we focus our attention on the R < 1 
case and we restrict to the case in which the system size 
is much shorter than the absorption length 

v gr (v gr Du 2 \ 

tabs = 7T~ = Ot a x . (61) 

27 \ C 7 e J 

The simplest quantity to measure in a static configura- 
tion is the time delay L/v gr : for sufficiently small group 
velocities, this allows the detection of very small differ- 
ences in L due, for example, to defects in the structure. 
Clearly this measurement requires a good temporal res- 
olution of the detector as well as some knowledge of the 
system parameters, in particular of v gr . 

The simplest example of dynamical modulation consists 
of a single-ramp modulation of il c . Only the small part 
of the pulse contained in the layer feels the modulation. 
As the layer is thin, the crossing time is often faster than 
the ramp time, which means that different parts of the 
pulse experience different portions of the ramp. In this 
case, the effect of the modulation on each given slice of 
the pulse can be estimated as 

1 n - n Vl, r . + Av nr (T) 

£/*(T = ^±->, (32) 

v 

where T is the time at which the slice exits the layer 
and the initial values v gr refer to the entrance of the 
slice in the EIT medium. If we approximate the ramp as 
linear, the variation in the electric field is then given by 
(Av/Vg r )(T/r), where Av is the amplitude of the group 
velocity ramp. 



As we have seen earlier in this article, a most interest- 
ing case consists of a double ramp, as illustrated in Fig. 
[6£a): the pulse is slowed down and then accelerated back 
to the initial group velocity. In the R < 1 case under 
consideration here, the part of the pulse which is mod- 
ulated during the slow-down ramp exits from the layer 
before the restoring ramp has begun. This latter ramp is 
then responsible for the creation of the peak in the trail- 
ing part of the pulse. The resulting shape is then very 
similar to the case shown in FigJH yet time reversed. 




8 10 12 14 



FIG. 6: Examples of wave-packet manipulation in several 
realizations of the EIT chain. All curves have been ob- 
tained using the model of Eq. (|29|l . Solid red lines show the 
modulated pulse, while dashed black lines result from static 
cases, (a) Single EIT layer [L = 200/^m, v gr = 10~ 7 c, 
T> = 6 x 1CP 7 (uip/kp)] with a Gaussian pulse (at = 10/xs, 
tabs ~ 500(7a;). Effect of a double ramp (as, e.g., in Fig[4|: 
v ~grl v gr = 10, t = 3.5 /is and storage time t s — 8 (is. (b) EIT 
double-layer structure [L = 30/im (each layer) , interlayer dis- 
tance AL = 3 x 10 7 /im, v = 5x 10~ 7 c, V = 3 x 10 -6 iv p /k'p] 
with a Gaussian pulse (a t = 1/is, £ a bs ~ 400 a x . Slow-down 
ramp in both layers: Av = —0.5 v gr , t = 50ns. (c) EIT chain 
with four layers [same as (b) but AL = 6 x 10 7 um] with Gaus- 
sian pulse [same as (b)]. Effect of a single slow-down ramp: 
Av — — 0.7 v gr , t — 50ns. 



B. Multilayer 

The single EIT layer is the basic building block for more 
complex geometries. As an example, the cases of single 
ramp in a double- and four-layer geometries are illus- 
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trated in Figs|6fb) andHJc). In particular, we restrict 
our attention to the case in which both the EIT and the 
vacuum layers are shorter than the pulse length. 

The simultaneous modulation of the dressing field on 
both layers allows the creation of several similar struc- 
tures on the same pulse, separated by a time depending 
on the distance between the layers. By choosing a fast- 
enough slow-down ramp, the peak that appears between 
the different layers can be shaped down to the absorption 
length. Its height can be enhanced by applying different 
ramps to the different layers, for example with opposite 
signs. The creation of such short and intense pulses can 
be of great interest in view of creating a strongly local- 
ized polariton, whose dynamics has been predicted to 
show peculiar features [201 ] . 



VII. CONCLUSION 

In conclusion, we have developed a Maxwell-Bloch for- 
malism to investigate the propagation of a light pulse 
through a multilayer structure consisting of alternating 
layers of a generic EIT medium and vacuum. The for- 
malism recovers the usual polariton dispersion in a ho- 
mogeneous EIT medium and is able to include reflection 
at interfaces as well as absorption processes. 

The effect of a dynamical modulation of the control field 
amplitude on the pulse shape is studied and specific at- 
tention is paid to the case when the modulation takes 
place while the pulse is overlapping an interface: as re- 



flection and losses remain negligible as long as the pulse 
fits in the EIT window, abrupt cuts on the wave-packet 
profile can be performed. 

In order to rapidly and efficiently simulate the propa- 
gation in geometry with strong velocity mismatch, an 
effective equation for the wave-packet propagation is de- 
veloped. The resulting modified continuity equation is 
shown to account for both propagation and absorption 
effects. A good agreement with the full Maxwell-Bloch 
calculations is found. 

Starting from this equation several manipulation schemes 
have been proposed and characterized using realistic pa- 
rameters for ultracold atomic clouds: contrary to stan- 
dard light-storage techniques, our proposal is able to 
take full advantage of the reduced size of typical atomic 
clouds. A most interesting possibility is the creation of 
highly localized peaks. 

Our calculations confirm that the wave-packet manipu- 
lation can be performed in a coherent and almost loss- 
less way, which opens interesting perspectives toward 
the quantum processing of single-photon wavefunctions 
for optical quantum computing applications. From a 
broader standpoint, the peculiar reflection properties at 
the interfaces of atomic media can be useful for polariton 
trapping and guiding as well as interface characterization. 
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